New Universality classes for generic higher character Lifshitz points 
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We describe new universality classes associated to generic higher character Lifshitz critical be- 
haviors for systems with arbitrary short range competing interactions. New renormalization- group 
arguments are proposed for anisotropic and isotropic systems. The usual m-axial Lifshitz universality 
classes are obtained as a simple limit of this arbitrary competing critical behavior. 



The universality classes characterizing the critical be- 
havior of competing systems are natural generalizations 
of those describing ordinary short-ranged critical sys- 
tems. In magnetic systems [1], the axial next- nearest- 
neighbor Ising (ANNNI) model [2] explains this gener- 
alized critical behavior, as antiferromagnetic exchange 
interactions between second neighbor spins are included 
along a single axis in a cubic lattice. The uniaxial Lifshitz 
point arises at the confluence of the disordered, a uni- 
formly ordered and a modulated phase. Second character 
Lifshitz points generalize the uniaxial one when the com- 
peting axes occur in arbitrary space directions. Physical 
applications of this behavior include magnetic modulated 
materials [3,4], high- Tc superconductors, liquid crystals, 
among others [5] . The theoretical description of this crit- 
icality has been put forth using different field-theoretic 
tools [6-8] and numerical Monte Carlo simulations for 
the ANNNI model [9]. 

The ANNNI model can be generalized by including 
further ferromagnetic couplings between third neighbors 
along the competing axis. The phase diagram is char- 
acterized by two parameters varying with the temper- 
ature, namely, the ratio of exchange couplings between 
the second and the first neighbors as well as the ratio 
of exchange interactions between the third and the first 
neighbors. The projection of the temperature axis in the 
plane of these parameters results in a region of intersec- 
tion where the different phases characterizing the system 
meet, known as a uniaxial Lifshitz point of third char- 
acter [10]. Including more alternate couplings along the 
competing direction, the number of parameters in the 
phase diagram increases. Analogously, the uniaxial Lif- 
shitz point of higher character is defined at the confluence 
of the distinct phases of the system. The Wi-fold Lifshitz 
point of character L arises when alternate couplings are 
permitted up to the Lth neighbors along directions 
and the ratios of the exchange interactions between dif- 
ferent neighbors are kept flxed at the associated Lifshitz 



temperature [11,12]. 

Lifshitz critical behaviors share at least two issues with 
a broad class of theories in other physical arenas. In some 
phenomenological models of quantum gravity, the Planck 
scale and the speed of light are two observer-independent 
scales which induce modifications of the usual energy- 
momentum dispersion relations [13]. A nonlinear real- 
ization of Lorentz invariance is maintained and can be 
generalized to include more than two relevant scales [14]. 
Lorentz invariance can also be preserved in addressing 
the dark matter problem in cosmology by making use of 
a scalar field with higher derivative terms, named "fc- 
essence" [15]. Alternatively, Planck scale effects may 
break Lorentz invariance in the ultraviolet regime due to 
the modification of the dispersion relations. Effective free 
quantum field theories with modified cubic kinetic terms 
have just been proposed [16] to test this idea. Thus, a 
renormalization group analysis with several independent 
length scales together with new perturbative techniques 
to solve interacting higher derivative field theories are 
highly desirable to elucidating new effects in a broader 
context. 

In this letter, we introduce the most general d- 
dimensional competing system resulting in a new static 
critical behavior by permitting only nearest neighbor 
interactions along (d — m2 — ... — rriL) directions, sec- 
ond neighbor competing interactions along m2 directions, 
etc., up to Lth neighbor competing interactions along ttil 
directions. The new universality classes for this system 
generalize the higher character universality classes in a 
nontrivial way. Since each competing direction defines 
distinct modulated phases, the generic higher character 
Lifshitz point arises at the confiuence among the disor- 
dered, a uniformly ordered and several modulated or- 
dered phases. The universality classes for this arbitrary 
competing exchange couplings Ising model, named here- 
after CECI model, depend on {N, d, m2, m^) where N 
is the number of components of the order parameter. The 
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particular isotropic universalities (rf = m„) arc obtained 
independently. We find multiscale scaling relations using 
a simple scaling hypothesis consistent with homogene- 
ity. Then, wc develop a new dimensional rcgularization 
procedure to calculate Feynman diagrams and critical ex- 
ponents for generic higher character Lifshitz points up to 
two-loop order. 

We start with the following bare Lagrangian density: 

n=2 

L ^ L—1 71—1 ^ 

+ 2 1 Vm„ M^ + Y.Y. 2 ' 

n=2 n=3 n'=2 

At the Lifshitz point, the fixed ratios among the ex- 
change couplings explained above translate into this field- 
theoretic version though the conditions (5o„ — r„„' = 0. 
All even momentum powers up to 2L become relevant in 
the free propagator, realizing a situation considered by 
Wilson in the 1970's [17]. 

The critical space dimension above which the clas- 
sical (mean field) behavior takes over the system is 
given by dc = 4 -f X]n=2 ^"^^^ ^w The small param- 
eter for a consistent perturbative expansion is = 

In momentum space, the m„ competing axes only pos- 
sess the term A;^^^ due to the Lifshitz condition. Thus, 
we can set cJn = 1 and perform a dimensional redefinition 



of the competing momenta defining arbitrary competing 
axes. Hence, if the canonical dimension (in mass units 
A) of the quadratic noncompeting momenta is given by 

[q] = A, then = A~. Therefore, the volume clement 
in momentum space for arbitrary Feynman integrals has 

canonical dimension A'' S»=2 n Let ^i, £,2t--7 

be the correlation lengths characterizing arbitrary spa- 
tial directions of the system. Close to the Lifshitz critical 
temperature these quantities are related among each, but 
can be considered as independent relevant length scales 
in the problem as follows. 

Consider L normalization conditions for the one- 
particle irreducible (IPI) vertex parts. The renormaliza- 
tion of the vertices at the Lifshitz temperature is consis- 
tent in the infrared limit if the theory is renormalized at 
nonvanishing external momenta. Let «;„ be the external 
momenta scale used to fix the renormalized theory along 
arbitrary space directions. For instance, for n = 1 the 
nonvanishing components of the external momenta are 
parallel to the noncompeting (ci— X^^=2 ?Tin)-dimensional 
subspace. For n > 1, the nonvanishing components of the 
external momenta occur only along the m„-dimensional 
competing subspace. Let Pj(„) be the external momenta 
associated to a generic vertex part. We can choose sev- 
eral symmetry points SPn along arbitrary directions de- 

2 

fined hj Pi(^n)-Pj{n) = ~" !)• The momentum scale 

of the two point fimction is fixed by p^^^ = = 1. 
The IPI renormalized vertex parts at the Lifshitz critical 
temperature are given in terms of the bare vertices as 
[8,18,19](n=l,...,L): 



R(n) \Piin) T ^ i(n) J 9n J l^n ) 



<P(n) («) (^\n) ' (Pi{n) , Qi(n) , A„ , A„) 



(2) 



-^Ar,o<^L,2r|°s^-'(Q(„), A„, A„)|q2 -^2)5 



where g„ are the renormalized coupling constants, .^^(,1) 
and 2'02(„) are the field and temperature normalization 
constants and A„ are the corresponding momentum cut- 
offs along the noncompeting (n = 1) and competing di- 
rections (1 < n < L). Notice that Qi(n) are the insertion 
external momenta along the different spatial directions 
for vertex parts including composite (^^(x) operators. 

Our dimensional redefinitions of the momenta compo- 
nents characterizing each type of competing axes implies 



that the renormalized and bare dimcnsionful coupling 
constants can be expressed in terms of dimensionless cou- 
plings as gn = Un{K^)^ and A„ = ■uo„(k^")"5^. Since 
we are going to analyze the situation for d < dc, we 
can suppress the cutoffs A„ in the argument of the bare 
vertices recalling, however, that the limit A„ — > 00 is im- 
plicit. The invariance of the bare vertex functions with 
the renormalizaton scales k„ implies the renormalization 
group (RG) equation 



dUn 



(3) 
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The functions /3n,70(n) and 702(„)(ii„) arc calculated at 
fixed bare coupling A„. In terms of dimensionless quan- 
tities, they read 



Pn = -neL[^ ) , 



dUn 



7<?i2(„)(Wn) = -Pn 



dlnZ, 



02 (n) 



dUn 



(4a) 
(4b) 
(4c) 



At the fixed points, under flows in the external momenta 

Pi{n)jQi{n) along arbitrary space directions, the vertex 

parts T^^(^^^\pPi(n), pQi(n),Ky i^n) havc a simple scaling 
property. As usual, define = 70(„)- The correlation 
length exponents are given hy = 2n — "y^2(^^y The 
determination of the scaling relations above and below 
the Lifshitz critical temperature follows closely that for 
the second character Lifshitz points [8,20]. One finds 



)!„ = 2 - n(rf - ^ ■ — mj)i'n, 

i=2 ^ 



(5a) 



(z-1) 



Pn = -^J^n{n{d ^ 
In = Vn{'^n-rin), 

n{d - Ef=2 ^^"ij) + 2n-r]n 



6n = 



nid - Ef=2 - 2n + r?„ 



rui) -2n + r]n), (5b) 
(5c) 
(5d) 



which imply the Widom 7„ = /3n('5n — 1) and Rushbrook 
an + 2Pn + 7n = 2 relations for arbitrary competing or 
noncompeting subspaces, since n = 1, ...,L. 

The isotropic cases are defined by d = m„, when all 
competing interactions not lying on the m„ competing 
axes are switched off. The upper critical dimension is 
dc = 4:n implying that the perturbative parameter is 
£4n = ^n — d. There is only one type of correlation length 
^4„. The canonical dimension of the momenta in mass 
units is redefined to be [k] = A^. The volume element 
of Feynman integrals in momentum space has canonical 
dimension [d'^'k] = A^. Following similar steps as done 
in the anisotropic cases, we obtain the relations: 



74n = Vin{2n - ri4n), 
a4n = 2 — dUin 

d + 2n- r]4n 



Sdn = 



•4n 



d - 2n + r]4n' 
^iy4n{d - 2n + r]4n) 



(6a) 
(6b) 

(6c) 
(6d) 



which imply the Widom 74„ = /34„((54„ — 1) and Rush- 
brook + 2P4n+^4n = 2 laws. Notc that they trivially 
reduce to the usual cj)^ theory for n = 1, after the iden- 
tification a4 = a, etc.. These relations are the same 



as those found by NicoU et al. [21]. with the additional 
scaling relation for /3„ obtained in the present work in a 
simple manner. 

Consider the exponents for the anisotropic cases ob- 
tained from the critical theory defined by (2) at nonzero 
external momenta. At this critical point the bare free 
critical propagator {t = 0) is given by G^^ = + 
Sn=2 ^n"' where p is a wave vector along the {d — TO2 — 
... — mi) noncompeting directions, and fc„ are the wave 
vectors along the m„ distinct type of competing axes. 
Let k[f^n) ^® external momenta along the m„ com- 
peting axes in the contribution of an arbitrary IPI ver- 
tex part, which at the symmetry points discussed above 
can be simply written as k'^^y If the loop momenta is 
A;(„), we can employ a generalized "orthogonal approx- 
imation" [8] to resolving the Feynman integrals up to 
two-loop level with arbitrary external momenta. Using 



the approximation (fc(n) + A:(„-))^" = {k'L)+^nY, 
redefine the loop momenta as K^^) = . In each inte- 
gration over the competing momenta we pick out only the 
contribution corresponding to a homogeneous function 
of the momenta. We absorb the geometric angular fac- 

tor [(n^^2 2k^™.r(f^))^(^_^.^^^^jr(2 - )] 

in a redefinition of the coupling constant after perform- 
ing each loop integral. Using (2), we find that the fixed 
point is indepent of the scaling direction considered and 
the critical exponents f„ and rjn can be calculated dia- 
gramatically at two- and three- loop order, respectively. 
The other exponents are obtained from the scaling laws 
derived above. They are given by 

1 _(iV+2)_ , 1 (jV-|-2)(Af2+23Af-f60) ,2 
_ _ 2 i{N+9.)''L T 8 (JV-|-8)3 x 

— — — 1 I ' a j 

n n 

n 2 ^ + 2 ,6(3iV+14) 1,. _ , 

'^" = "^^ = r^(iVT8F[^ + ^"( (Ar + 8)^ -4^^' 



we can 



+ 2(7V + 8)' 

(4-^) _ (iV + 2)(iV^ + 30iV + 56) ^ 
2{N + 9,f^ 4(7V + 8)3 ^ > 

B - 1 3 (jV + 2)(2iV+l) „ 

^""2 2(7V + 8)''^ 2(7V + 8)3 

Sn-3 + e,+ 2(7V + 8)2 ^^^^ 

Notice that the anisotropic exponents 7n,a„,/3„ and 5„ 
easily reproduce those of the second character Lifshitz 
behavior when = ... = rriL = with m2 = m ^ 0, 
which turn into the Ising-like exponents for m2 = m = 0. 
The generalized universality classes exposed here and the 
properties i/„ = ^ and r]^ = rir}i, when m2 = ... = 
rriL-i = with ^ 0, imply that our anisotropic scal- 
ing relations are equivalent to those found in Ref. [12], 
except for the scaling relation of the magnetization ex- 
ponent which was not obtained there. Strong anisotropic 



2) (Ar + 2)(A^2 + 22A^ + 52) ^ . 

+ 4(7v + 8)3 



Oin = 



3 



scaling [22] is exact at the loop order considered here in 
the most general anisotropic situation, and might be ex- 
pected to hold at arbitrary loop levels. 

Using the same approximation scheme to perform 
Feynman integrals as outlined above, the isotropic ex- 

I 
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OLAn 
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The isotropic exponents q:4„, /?4„, 74^ and Sin depend ex- 
plicitly on the number of neighbors coupled through com- 
peting exchange interactions n and labels the isotropic 
universality classes. These critical exponents reduce to 
those from the Ising-likc {(f)^) miivcrsality class for n = \. 
For n > 2, these exponents can never be obtained from 
those of the anisotropic cases described above, since the 
latter do not depend explicitly on n. This confirms that 
the isotropic universality classes are independent from 
the anisotropic universality classes for Lifshitz points of 
generic higher character. 

In summary, the CECI model proposed here presents 
the most general sort of second order phase transition for 
arbitrary types of competition and realizes new univer- 
sality classes for generic higher character Lifshitz critical 
behaviors. The scaling relations are obtained through a 
generalized scaling hypothesis which defines independent 
renormalization group tranformations in each competing 
(or noncompeting) subspace. The critical exponents are 
calculated up to two-loop level using a new analytical 
dimensional regularization method to compute Feynman 
integrals with arbitrary even momentum powers in the 
propagator, which may be fruitful in quantum field the- 
ory applications. Universal quantities as amplitude ratios 
[24,25] can be pursued in this framework. We relegate 
the details of the results described here for a subsequent 
publication [26]. 

It is important to mention that a new anisotropic be- 
havior emerges when d = m2 + ... + "/l. The noncom- 
peting directions which tend to favor the existence of the 
ferromagnetic phase, are absent in this new situation. 
This particular behavior deserves a complete and careful 



ponents are even simpler to obtain. The angular factor 
to be absorbed in the loop integrals is simply the area of 
the m„-dimensional sphere Smr,- The critical exponents 
are given by [23]: 

(8a) 
(8b) 
(8c) 
(8d) 
(8e) 
(8f) 

study for future investigations. The analysis of the most 
general crossover effects as well as the tricritical behavior 
of the model are also worthwhile. 
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